Critical spin liquid at 1/3 magnetization in a spin- 1/2 triangular antiferromagnet 
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Although magnetically ordered at low temperatures, the spin- 1/2 triangular antiferromagnet CS2CUCI4 ex- 
hibits remarkable spin dynamics that strongly suggest proximity to a spin liquid phase. Here we ask whether a 
proximate spin liquid may also occur in an applied magnetic field, leaving a similar imprint on the dynamical 
spin correlations of this material. Specifically, we explore a spatially anisotropic Heisenberg spin- 1/2 triangu- 
lar antiferromagnet at 1/3 magnetization from a dual vortex perspective, and indeed find a new "critical" spin 
liquid phase described by quantum electrodynamics in (2+l)-dimensions with an emergent SU(6) symmetry. A 
number of nontrivial predictions follow for the dynamical spin structure factor in this "algebraic vortex liquid" 
phase, which can be tested via inelastic neutron scattering. We also discuss how well-studied "up-up-down" 
magnetization plateaus can be captured within our approach, and further predict the existence of a stable gap- 
less solid phase in a weakly ordered up-up-down state. Finally, we predict several anomalous "roton" minima 
in the excitation spectrum in the regime of lattice anisotropy where the canted Neel state appears. 

PACS numbers: 



The quest for an unambiguous experimental realization of 
a quantum spin liquid remains a central pursuit in condensed 
matter physics, despite a long history dating back to An- 
derson's suggestion that the nearest-neighbor Heisenberg tri- 
angular antiferromagnet may realize a "resonating valence 
bond" ground state. Jj]] Quite generally, the theoretical search 
for models realizing such exotic quantum ground states has fo- 
cused primarily on frustrated magnets in zero magnetic field. 
The central goal of this paper is to take a first step at analyz- 
ing the situation when a finite magnetic field is present and 
ask to what extent spin liquids may occur in this broader set- 
ting. Naively speaking, this may seem somewhat misguided, 
since sufficiently strong magnetic fields quench quantum fluc- 
tuations entirely and lead to a simple ferromagnetic ally or- 
dered ground state. However, at moderate-strength fields it is 
conceivable that the presence of numerous competing phases 
arising from the geometric frustration may lead rather to an 
enhanced role of quantum fluctuations by the field. 

The spin- 1/2 triangular antiferromagnet provides a sim- 
ple and experimentally relevant test case for these ideas. 
One noteworthy example is the spatially anisotropic material 
Cs2CuCl4,||2l |3l |4[], which has garnered much attention as a 
promising experimental realization of a spin liquid at zero 
magnetic field. Though CS2CUCI4 exhibits long-range spiral 
order at the lowest temperatures, inelastic neutron scattering 
experiments reveal broad regions of continuum scattering at 
intermediate energies throughout the Brillouin zone.j4|] This 
anomalous scattering coexists with well-defined spin-waves 
in the ordered phase and, significantly, also persists at tem- 
peratures above the Neel temperature where the spin waves 
are absent. The origin of this continuum scattering presents 
one of the foremost challenges for understanding the behav- 
ior of this interesting material. While two groups find that 
nonlinear spin-wave theory can account for much of the ob- 
served weight in the ordered phase, spin-liquid physics 
has been widely invoked as a possible explanation for the 
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Low-temperature magnetic order develops in the pres- 
ence of a magnetic field as well, leading to a rich phase 
diagram. [15] Motivated in part by the observed zero-field phe- 
nomenology, we will explore the following question here. Can 
spin liquid phases appear at non-zero magnetic field which in- 
fluence the dynamics of CS2CUCI4 at intermediate energies, 
just as appears to be the case in the absence of a field? Several 
experimental features] 15] are worth noting that make this sce- 
nario plausible. First, a broad temperature range characterized 
by short-range order persists up to sizable fields of around 6 T. 
Second, the ordering temperature initially decreases with the 
magnetic field strength, thereby broadening this short-range 
order regime. Finally, there is a stark contrast in the experi- 
mentally determined phase diagrams for moderate fields ap- 
plied along the b and c axes in the plane of the triangular lay- 
ers, implying a high-sensitivity of ordering to small perturba- 
tions. These features together strongly point to the presence of 
many nearly degenerate states and a corresponding enhance- 
ment of quantum fluctuations by intermediate-strength fields. 

To explore possible field-induced spin liquid phases rele- 
vant for CS2CUCI4, we study an anisotropic Heisenberg spin- 
1/2 triangular antiferromagnet tuned to 1/3 magnetization us- 
ing a well-studied duality mapping! 16]. The present paper 
significantly extends an earlier study of an (easy-plane) spin 
model in the absence of a magnetic field. 1 1 ill At zero-field, 
it was argued that reformulating the spin model in terms of 
fermionized vortices leads naturally to a spin liquid of the 
"critical" (versus topological) variety. This "algebraic vor- 
tex liquid" (AVL) is a promising candidate for explaining the 
zero-field CS2CUCI4 phenomenology. As we describe below, 
generalizing to the case of 1/3 magnetization again leads nat- 
urally to a "critical" AVL phase, which like its zero-field pre- 
decessor supports gapless spin excitations and power-law spin 
correlations. The new AVL characterized here is described by 
quantum electrodynamics in (2+l)-dimensions (QED3) with 
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an emergent SU(6) symmetry, which has important implica- 
tions for the spin dynamics. Specifically, it follows that both 
the components of the dynamical spin structure factor along 
the field and perpendicular to it exhibit enhanced universal 
correlations with identical power-laws at various momenta 
in the Brillouin zone (i.e., the filled circles in Fig. |2j. Such 
anomalous scattering would be quite interesting to search for 
in CS2CUCI4 via inelastic neutron scattering. On a technical 
level, we emphasize two appealing features of this AVL. First, 
the field explicitly breaks the SU(2) spin symmetry down to 
U(l), so that our dual description is valid even in the absence 
of a Dzyaloshinskii-Moriya interaction, thereby avoiding sub- 
tleties encountered at zero field. Second, the QED3 theory 
describing this AVL is "larger- N" than in the zero-field case, 
and thus even more likely to exist as a stable phase. 

We also discuss the so-called "up-up-down" (UUD) 1/3 
magnetization plateaus from our dual perspective. Such 
UUD states have been well-studied for the isotropic trian- 
gular antiferromagnet using spin-wave theory ! 1711 and exact 
diagonalization| 18], and have recently been observed in the 
anisotropic material Cs2CuBi4[19]. We show how the 1/3 
magnetization plateau can be captured within our approach 
when the sublattice magnetization is near full polarization, 
and further predict the presence of a "critical" solid when the 
sublattice magnetization is weak. Signatures of this gapless 
solid phase may be accessible in exact diagonalization and/or 
series expansion studies by examining the excitation spectrum 
in an XXZ model with spin-space anisotropy tuned near the 
transition to the UUD plateau. 

Finally, we discuss the range of lattice anisotropy where the 
ground state is the canted square-lattice Neel phase. In this 
"frustrated square lattice" limit, we predict several anomalous 
minima in the excitation spectrum at the momenta indicated 
by open circles in Fig. [2] which are due to vortex-antivortex 
"roton" excitations. These excitations can be difficult to cap- 
ture within spin-wave theory but appear naturally in our dual 
framework. Such features were also predicted at zero-field, 
in agreement with earlier series expansions. Il20ll and would be 
interesting to probe here as well to (perhaps) further substan- 
tiate the vortex interpretation of those results. 

Turning to the details, the Hamiltonian we consider is 



H = 



>J rr' St- * St" 
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where the magnetic field h lies in the triangular planes and the 
anisotropic exchanges are as shown in Fig.[T] Throughout we 
assume that the field is tuned so that the system is at 1/3 mag- 
netization. Furthermore, we will ignore the small interplane 
and Dzyaloshinskii-Moriya couplings present in CS2C UCI4 , 
which is appropriate for the field orientations of interest 02lll . 
Since the Hamiltonian has U(l) spin symmetry, Eq. (1) can be 
mapped onto a quantum rotor model and dualized using the 
standard transformation [16]. Readers interested in the details 
of this transformation and the analysis to follow are referred to 
Refs. Il22ll and 01211 . where the approach we employ has been 
well-developed in similar settings. Here we will apply this 



technique to the new physical situation of a 1/3-magnetized 
triangular antiferromagnet, highlighting the nontrivial physics 
it enables us to access but dispensing with unnecessary for- 
malism developed elsewhere. 

The quantum rotor mapping is implemented by replacing 
5+ — > e liPl and S z — * n r — 1/2, where n T is an integer- valued 
boson number and cp T is the conjugate phase. While not exact, 
such a transformation is expected to be inconsequential for 
describing the universal physics that is our focus. The rotor 
Hamiltonian can then be expressed as 

H' = ^J rr >cos{(p r -(p r ,) + U^2{n r -l/3) 2 

(rr') r 

+ Jrr'K-l/3)(rar' -1/3), (2) 

(rr') 

where the U term above enforces energetically the constraint 
of having either or 1 boson per site as appropriate for mod- 
eling a spin- 1/2 system. In the quantum rotor language, the 
condition of 1/3 magnetization translates into having on aver- 
age one boson for every three sites, which is manifest in the 
above Hamiltonian. 

The duality mappin g ap plied to Eq. (fJJ proceeds in an iden- 
tical fashion as in Refs.l22landll2[ the only difference between 
these references and the present system is that the bosons are 
now at a different mean filling. In the dual picture, one equiv- 
alently reformulates the rotor model in terms of quantum me- 
chanical, bosonic vortex degrees of freedom, which are topo- 
logical defects in which the spins wind around triangular pla- 
quettes as shown in Fig. Q] These vortices are mobile, point- 
like particles that hop on the dual honeycomb lattice (see Fig. 
03 in a background of a fluctuating gauge field a XX ' whose 
flux encodes the boson number (or equivalently the S z com- 
ponent of spin, along the field), n z ~ (A x a) r /(27r). Thus, 
the magnetic field manifests itself as a nontrivial background 
flux "felt" by the vortices, which at 1/3 magnetization is a 
commensurate 2tt/3 flux per dual hexagon on average. This 
background flux is where the present study departs from the 
zero-field analysis of Refs. IjjJ, [l2j, where the vortices see 7r 
flux, and is responsible for the new physics we obtain here. 
The vortices interact via a logarithmic repulsion mediated by 
the gauge field, and importantly are at half-filling due to the 
underlying frustration in the original spin model. In terms of 
a vortex number operator N x and its conjugate phase 6 X , the 
dual Hamiltonian takes the form 

Hduai = — ^ t xx > cos(9 x — 9 x i — a xx i) 

<xx'> 

+ ^(Wx - 1/2)V XX ,(N X , - 1/2) + H a . (3) 

xx' 

The first term allows vortices to hop across nearest-neighbor 
honeycomb sites. The hopping amplitudes are generally 
anisotropic since vortices hop more easily across weak spin 
links; thus we take t'/t ~ J' / J (see Fig. [TJ. The second 
term encodes the vortex repulsion, while the last describes the 
gauge field dynamics. 
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FIG. 1 : (Color online) Triangular lattice and the dual honeycomb on 
which vortices reside. Spins shown illustrate a vortex. 



The large vortex density poses a significant challenge for 
analyzing the dual theory as it stands. The strong vortex inter- 
actions, however, actually make the problem more tractable, 
as they strongly suppress vortex density fluctuations and lead 
to an incompressible vortex fluid. Vortex exchange statistics 
is therefore of secondary importance, and one can proceed by 
utilizing a formally exact mapping to convert the bosonic vor- 
tices into fermions bound to 2ir flux tubes. This statistical flux 
does not alter the average flux seen by the vortices, since it 
averages to 2ir per hexagon (which is equivalent to zero flux). 
Furthermore, as argued in detail in Ref.[ljfhe incompressibil- 
ity renders the flux attachment irrelevant for describing low- 
energy physics of the vortex fluid. This statistical transmuta- 
tion at low energies from bosonic to fermionic vortices is at 
the heart of our approach. 

Physically, working with fermions is advantageous because 
the Pauli principle allows one to first focus on the vortex ki- 
netic energy while still maintaining a good interaction en- 
ergy. Thus, we initially consider a mean-field state where 
we "smear" the 2ir/3 background flux (due to the 1/3 mag- 
netization) uniformly across the lattice, though we will relax 
this assumption later. Fluctuations around this flux-smeared 
state can be systematically controlled as discussed below. 
Within this mean-field, one can work out the vortex band 
structure, which captures the important intermediate length- 
scale physics. In particular, we find that for t'/t > 2 1 / 3 the 
vortices are gapped, while for t'/t < 2 1 / 3 they form a "crit- 
ical" state with six gapless Dirac points. The former gapped 
state corresponds to the canted "square lattice" Neel phase ex- 
pected when J' is dominant, while the latter is the mean-field 
description of the AVL which will be our main focus. 

We discuss the gapless regime first. Expanding around 
the Dirac points and including gauge fluctuations around the 
flux-smeared state, one obtains a low-energy effective theory 
describing six flavors of two-component Dirac fermions ip a , 
a = 1, . . . , 6, which are minimally coupled to a U(l) gauge 
field a M that mediates the vortex repulsion. The low-energy 
effective theory so obtained is identical to QED3, 

£QED3 =^ a {9~ + -^{t^\dva\) 2 + £ 4 f. (4) 

(Inclusion of the statistical 2tt flux tubes merely leads to inter- 



actions that are irrelevant by power-counting in the above, as 
claimed earlier.) The first term encodes the linearly dispersing 
kinetic energy for our six flavors of vortices. The second is the 
usual Maxwell term, while represents symmetry-allowed 
four-fermion interactions. Of central importance is whether 
the "critical" nature of the mean-field state survives in the full 
interacting theory above. A partial answer to this question 
can be obtained by asking whether any fermion mass terms, 
which would drive some type of ordering in the original spin 
model and lead to a gap in the vortex spectrum, are allowed 
by symmetry. Despite the loss of time-reversal symmetry by 
the magnetic field, the answer is "no" — the remaining symme- 
tries of the original spin model are still sufficient to preclude 
all possible mass terms in Eq. ©. This is not the full story, 
however, since the four-fermion interactions above, if relevant 
in the renormalization group sense, could still potentially lead 
to ordering via spontaneous mass generation. Hence to pro- 
ceed we must assess the role of these terms in the theory. 

In the limit of a large number N of fermion flavors, all 
such four-fermion interactions are indeed known to be irrel- 
evant, so that QED3 realizes a nontrivial stable critical phase. 
(See QED3 refs. in 11211 .) While the critical value of N above 
which this holds is uncertain, calculations to leading order in 
1 /N suggest that N — 6 relevant here is large enough. Hence, 
we proceed with the assumption that Eq. indeed describes 
a stable critical phase for our vortices, with an emergent SU(6) 
symmetry due to the presumed irrelevance of C±{ . In terms of 
the original spin model, this can be qualitatively understood as 
follows. The presence of numerous gapless Dirac points im- 
plies that there are many competing orders in the spin model. 
With sufficiently many competing orders {i.e., at large enough 
"N"), quantum fluctuations can be so strong as to disorder 
the system even at zero temperature. The resulting critical 
phase is precisely the AVL, which respects all symmetries of 
the original spin model and supports gapless vortex excita- 
tions and, in turn, gapless spin excitations as we now discuss. 

The key experimental prediction for this phase lies in the 
behavior of the dynamical spin structure factor, since this can 
be directly probed with inelastic neutron scattering. We first 
discuss the spin correlations of S^, transverse to the field. 
Recalling that S z + 1/2 ~ (Ax a)/(2vr), since S + adds 
S z = 1 it follows that the corresponding dual operators are 
"monopoles" which add 2tt gauge flux. The added flux gives 
rise to six additional vortex zero-modes, one for each fermion 
flavor, and half of these must be filled to produce a physical 
state. Thus there are 20 leading monopole excitations, which 
can be shown to carry the momenta IL, displayed in the left 
side of Fig.|2][12] Such monopoles exhibit nontrivial power- 
law correlations, each with identical scaling dimension due to 
the emergent SU(6) symmetry. Consequently, near each IT,, 
the transverse spin structure factor scales as 
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{cu 2 - q 2)i-u + _/a- 



(5) 



Using the leading large- TV result. B23I1 the anomalous dimen- 
sion is ?7 + _ w 0.54N - 1 « 2.2 for each IT,-. 
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FIG. 2: Filled circles indicate momenta at which the components of 
the spin structure perpendicular to the field (S + ~) and along the field 
(S zz ) exhibit enhanced universal correlations with the same power- 
law decay in the AVL. In the canted Neel phase, anomalous "roton" 
minima in the excitation spectrum are predicted at the momenta de- 
noted by open circles. 



Consider now the correlations of S z , along the field. Near 
zero-momentum S z appears as the dual gauge flux Ax a. But 
the vortex band structure also allows naturally for "particle- 
hole" excitations which correspond to vortex currents, and 
such currents generate modulated gauge flux that contributes 
to S z at other wave vectors. These appear in the continuum 
as fermion bilinears, which provide the leading S z spin cor- 
relations at the 18 momenta denoted by filled circles in the 
right side of Fig. [2] Near each of these momenta the struc- 
ture factor S zz scales as in Eq. ©, but with a larger anoma- 
lous dimension r/ zz w 2.3 estimated from the leading 1/JV 
result. [24] Note that the momenta denoted by filled circles in 
Fig. |2] shift along k x as the lattice anisotropy changes, lining 
up along the dashed lines in the 1-D limit and residing sym- 
metrically around the hexagons shown in the isotropic limit. 

We also mention that both r;_| and r\ zz are estimated to be 

larger than 2, in which case a cusp rather than a divergence oc- 
curs in the structure factor. Enhanced scattering should never- 
theless be observable near the above momenta, and moreover, 
in CS2CUCI4 such divergences would be cut off at the lowest 
energies due to the onset of magnetic order. 

We now shift our attention to 1/3 magnetization plateaus, 
focusing for simplicity on the isotropic limit J' = J. Return- 
ing to the "flux-smeared" mean-field level, we now weakly 
modulate the gauge flux away from 27r/3 such that UUD or- 
der is assumed at the outset. Introducing this flux modulation, 
surprisingly, merely shifts the locations of the Dirac nodes. 
Furthermore, the remaining symmetries of the spin model are 
still sufficient to protect the gaplessness of this state against 
small perturbations. That is, the critical nature of the AVL 
is preserved upon introducing weak UUD order, implying the 
interesting possibility of having a stable gapless solid phase. 

Where, then, is the gapped magnetization plateau! Such 
a spin state is described by a v = 1 quantum Hall state 
for the fermionic vortices. To see this, note that the 
first quantized wave function for the bosonic vortices is 
vp bos = JT e~ ie(Xl ~ Xj) *fcrm, where 0(x) denotes the an- 
gle formed by the vector x with respect to a fixed axis and 



'I' form is the fermionic vortex wave function. With 'J/ferm 
the usual v = 1 wave function, it has been shown|25[ that 
^bos exhibits (quasi) off-diagonal long-range order. Hence 
the bosonic vortices form a "superfluid", which corresponds 
to an "insulating" state for the original spin model. Indeed, 
due to the accompanying dual "Meissner effect" the gauge 
field a XX ' picks up a Higgs mass, and there are thus no gap- 
less excitations in this phase. This is the desired UUD plateau 
as claimed. Such a v = 1 quantum Hall state can be driven 
by adding an imaginary second neighbor vortex hopping with 
strength ti larger than a critical value t^c, which decreases as 
the sublattice magnetization increases towards full polariza- 
tion. This can be explicitly demonstrated by computing the 
Chern numbers for the occupied bands. [26] Second neighbor 
hoppings with t 2 < t 2c lead to symmetry breaking, and are 
thus precluded in accordance with the above discussion. 

Finally, let us discuss the "frustrated square lattice" limit 
where J' is dominant and the vortices are gapped. Again, 
this regime corresponds to the expected canted Neel phase. 
Consider the correlations of S z at momentum q. Aside from 
spin-waves, the spin structure factor S zz receives contribu- 
tions from vortex-antivortex "roton" excitations. In the spin 
language, these correspond to vortex currents generating mod- 
ulated S z as discussed above in the AVL. Such excitations are 
analogous to Feynman's rotons in He-4, and likewise should 
appear as minima in the structure factor. The energy required 
to create a roton with momentum q is simply given by the 
minimum energy required to promote a fermionic vortex in an 
occupied band with momentum k to an unoccupied band with 
momentum k — q. lll2ll As J increases, enhancing frustration, 
the vortex band gap shrinks leading to a sharp reduction in the 
minimum roton energy at the commensurate wave vectors de- 
noted by open circles in Fig. [2] When the gap closes signaling 
the destruction of the canted Neel order, the roton energy be- 
comes gapless at these momenta. (For larger J one enters the 
AVL phase, and the additional momenta in the right side of 
Fig. |2] denoted by filled circles then branch out from these ro- 
ton minima.) The presence of these low-energy rotons should 
lead to dramatic deviations from linear spin-wave theory. 

To conclude, we have provided the first concrete theo- 
retical proposal for a spin-liquid which may influence the 
intermediate-energy dynamics of Cs 2 CuCl4 in a magnetic 
field. The prospect of observing the spin-liquid physics de- 
scribed here is an exciting one, and we hope experiments in 
this direction will be pursued. Our nontrivial predictions for 
the dynamic spin structure factor can be tested with inelastic 
neutron scattering by measuring the lower-edge of the contin- 
uum scattering at the momenta specified in Fig. |2] Polarized 
neutrons in particular would provide a useful probe for the 
markedly different correlations identified parallel and trans- 
verse to the field. We also identified a new stable gapless 
phase with weak UUD order, which would be interesting to 
search for via exact diagonalization and series expansions. A 
renewed look at the excitation spectrum in the UUD plateau 
as one adds easy-plane anisotropy to suppress the solid order 
may prove fruitful. Series expansion studies to search for the 
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predicted rotons in the "frustrated square lattice" limit with J' 
dominant would also be interesting. More generally, the new 
spin-liquid presented here suggests that it may be worthwhile 
to widen the search for such exotic phases in other frustrated 
systems by incorporating a finite magnetic field. 
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